Coexistence of d I .2_ !/ 2-wave superconductivity and antiferromagnetism 

induced by a staggered field 
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The two-dimensional t-J model in a staggered field is studied by exact diagonalization of small 
clusters. For the low-hole-density region and a realistic value of J/t, it is found that the presence of 
a staggered field strengthens the attraction between two holes. With increasing field, the d x 2_ y 2- 
wave superconducting correlations are enhanced while the extended-s-wave ones hardly change. 
This implies that coexistence of the d I ,2_ !/ 2-wave superconducting order and the commensurate 
antiferromagnetic order occurs in a staggered field. 

PACS numbers: 74.20.-z, 74.20.Mn, 71.27. +a 



I. INTRODUCTION 

In high-T c cuprates, the antiferromagnetically ordered 
phase and the superconducting phase were separately ob- 
served in the plane of doping concentration and tempera- 
ture. 1 ' 2 Thus antiferromagnetism and superconductivity 
may be thought as competing with each other. However, 
some of the recent experiments suggest a possibility of 
their coexistence, although it is still controversial. 3-8 

Elastic neutron scattering experiments in YJi^Cu^Oy 
with y = 6.5 and 6.6 show that magnetic intensity 
emerges near room temperature at the momentum (jr, it) 
in units of the reciprocal lattice parameter. 3,4,7 The in- 
tensity increases continuously with decreasing tempera- 
ture. Remarkably, an upturn of the intensity is observed 
at a superconducting transition temperature. This sug- 
gests coexistence of superconductivity and the commen- 
surate antiferromagnetic (AF) order. Moreover, nu- 
clear quadrupole resonance measurements have revealed 
the presence of magnetic moments in the superconduct- 
ing state on Hgo.8Cuo.2Ba2Ca2Cu3C , 8+ ( 5- 8 The observed 
magnetic moments in these materials are often regarded 
as a consequence of the formation of the d-density wave 
order. 9 However, this interpretation leaves a difficulty 
since the d-density wave order suppresses superconduc- 
tivity. 10 Thus we should also consider an alternative sce- 
nario that the magnetic moments are due to the ordered 
Cu spins. 

So far, the possible coexistence of antiferromagnetism 
and d-wave superconductivity in the two-dimensional 
(2D) t-J model has been discussed at low hole doping, 
in a variational approach 11 and quantum Monte Carlo 
calculations. 12 The possibility of the coexistence 13 poses 
a fundamental question on interplay between antiferro- 
magnetism and d-wave superconductivity. Before an- 
swering whether the coexistence actually takes place in 
cuprate superconductors, we would like to clarify whether 
those two orders can coexist in strongly correlated elec- 
tron systems. To clarify the matter, it would be useful to 
study the hole pairing and superconductivity in a stag- 
gered field, which forces the system to have the AF or- 
der. Indeed, for the ID t-J model in a staggered field, 



the superconducting correlation was found to be the most 
dominant for J/t ~ 0.4. 14,15 

In this paper, we investigate the 2D t-J model on a 
square lattice in a staggered field coupled to electron 
spins. While the staggered field is introduced here as an 
artificial parameter to induce the AF order, it may arise 
naturally if three-dimensional inter-plane interactions are 
treated in a mean field theory. We employ exact diag- 
onalization for the 4x4, x VT8, y/20 x y/20, and 
V26 x V26 clusters with periodic boundary conditions. 
Our results demonstrate that a staggered field actually 
enhances the pairing of two holes and the d^^^-wave 
superconductivity. 



II. MODEL 

We consider the following Hamiltonian given by 



(i,J) 



(1) 
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where is the nearest neighbors. The constrained 

fermion operator cj CT is given by = c^(l - 
which means that double occupancy at each site is ex- 
cluded. The last two terms are due to the presence of a 
staggered field whose magnitude is denoted by h; A and 
B represent the two sublattices on a square lattice. We 
refer to this model as the t-J-h model. In this work we 
fix J/( = 0.4 which is considered as a realistic value, and 
vary h/t as a parameter. 



III. NUMERICAL RESULTS 

We first discuss the hole correlation function given 
by C hole (f) = (l/N)J2t{n h (i)n h (t + r)). Here N is 

the number of lattice sites, ?ih(i) = 1 — n*, and (• • •) 
denotes the expectation value in the zero-momentum 
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FIG. 1: Left: Equal-time hole correlations as a function of dis- 
tance in the 2D t-J-h model with N = 18, J/t = 0.4, and var- 
ious values of h/t. Crosses, circles, triangles, and squares are 
the data for h/t = 0.0, 0.4, 1.0, and 2.0, respectively. Right: 
Root-mean-square separation of the hole pair as a function of 
h/t. J/t = 0.4. 



ground state. In the left panel of Fig. 1 we show the 
distance dependence of CWcM with r = \ f] for the hole 
densities N h /N = 2/18 ~ 0.111 and 4/18 ~ 0.222. For 
two holes and h = 0, the most dominant correlations are 
at r = \/2, namely, when the holes stay at the next- 
nearest neighbors. 16 ' 17 As h/t increases, correlations at 
the nearest neighbors (r — 1) become stronger than those 
at r = V2, and contribution at longer distances is sup- 
pressed. For four holes, correlations at r = 1 are en- 
hanced while ones at the largest distance hardly change. 
This means that the presence of a staggered field makes 
the interaction between two holes attractive but the hole- 
pairs are well separated. 

In order to analyze the obtained data in the two- 
hole case, we calculate the root-mean-squar e sep aration 
of the hole pair 18,19 defined as r rms = \J (r 2 ) where 

(r 2 ) = E^Pl 2 Choie(0/E^3)Choie(rO. Here \f'\ 
takes the shortest distance between two holes on the lat- 
tice with periodic boundary conditions. The right panel 
of Fig. 1 shows the /i/i-dependence of r rms . It is clear that 
the separation of the hole pair becomes smaller with in- 
creasing staggered field. For larger h/t, the value of r rms 
is less affected by finite-size effects. The present result 
suggests that the staggered field does help binding of two 
holes, and that the hole binding survives in the thermo- 
dynamic limit at least for a sufficiently large h/t. 



A tendency of hole binding can be obtained also from 
the binding energy, which is given by 20 

E B = E (N h = 2) + E Q (N h = 0) - 2E Q (N h = 1). (2) 

Here E (Nh) denotes the ground-state energy with N-^ 
holes in N sites. A negative value of indicates the 
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FIG. 2: (a) Binding energy as a function of h/t in the 2D t-J- 
h model with J/t = 0.4. (b) Size dependence of the binding 
energy. Crosses, pluses, circles, and triangles are the data for 
h/t = 0.0, 0.2, 0.4, and 1.0, respectively. J/t = 0.4. 



presence of hole binding. In Fig. 2(a) we show the depen- 
dence of the binding energy on the staggered field h. The 
binding energy is negative in the whole range of h/t and 
has a peak at h/t ~ 0.8, 1.2, 1.2, and 0.8 for N = 16, 18, 
20, and 26, respectively. The obtained results apparently 
imply that hole pairing is suppressed by a small staggered 
field. However, being an intensive quantity, the binding 
energy is severely affected by finite-size effects. 18 Figure 
2(b) shows the size dependence of the binding energy for 
various values of h/t. In fact, although at h = (and 
J/t = 0.4) the binding energy for N < 26 takes negative 
values, an extrapolation from the data rather indicates 
the absence of hole binding in the thermodynamic limit. 
This has been already discussed in earlier studies 21,22 for 
N < 32. On the other hand, for larger h/t, we find 
the size dependence to be substantially smaller. It seems 
that the binding energy remains negative with increasing 
iV in the presence of a staggered field. Therefore there 
is a possibility that two holes in the bulk limit tend to 
be bound even by a small field. In particular, for a large 
h/t (h/t <; 1.0) the weak size dependence of the negative 
binding energy strongly suggests the hole pairing in the 
bulk limit. 



The pairing of holes is also consistent with the 
enhanced superconducting correlation discussed below. 
We calculate the equal-time superconducting corre- 
lations given by C a (r) = (l/N) £ ? <At(*)A Q (t + 
rj). 23,24 The singlet pairing operator A a (i) is defined as 

A Q (t) = (1/V2) E?/«(e) ( c ? T c i+r,i ~ % c U?a) > where e " 
is (±1,0) and (0, ±1). For the extended-s-wave pairing 
symmetry (a = s), wc put / s (e) = +1 at all e. For 
the d x 2_ y 2 symmetry (a = d), we put /d(e) = +1 at 
e = (±1,0) and f d (e) = -1 at e = (0,±1). Figure 3 
shows the distance dependence of Cd (r) and C s (r) for the 
hole densities 0.111 and 0.222 in 18 sites. For n h ^ 0.2, 
with increasing field, the d x 2_ y 2--w&ve superconducting 
correlations are enhanced at all distances with r > 1. In 
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FIG. 3: Equal-time superconducting correlations as a function 
of distance in the 2D t-J-h model with N = 18, J/t = 0.4, 
and h/t = 0.0, 0.4, 1.0, and 2.0. The symbols are the same as 
in the left panel of Fig. 1. 



contrast, the extended-s-wave ones hardly change, espe- 
cially at long distances. This implies that the presence 
of a staggered field helps the d^.^-wave superconduc- 
tivity. 



Calculation of the pair spectral function 25 ' 26 should 
provide another evidence for the d a; 2_ y 2-wave pairing en- 
hanced by a staggered field. The pair spectral function 
is given by 

P a (uj) = ^2\^n(N h = 2)\A t : t \^ o (N h ^0))\ 2 

n 

x5{lo - E n (N h = 2) + E (N h = 0) + /x),(3) 



where A** = E ? A a (i)/V]V, \i = E Q {N h = 2) - 
Eo(Nh — 0), and \^ n (Nh)} denotes an eigenstate with 
energy E n (N\ 1 ) in the A^-hole system. The left panel 
of Fig. 4 shows the w-dependence of Pd(w) for various 
values of h/t. The overall feature approximated by a 
Lorentzian is insensitive to the system size at each h/t. 
The peak at 10 — (i.e., the coherent peak) grows with 
increasing h/t, which means that the pairing becomes 
strong. The contribution for w > 0, which seems to be 
a continuum spectrum, is relatively suppressed, but a 
peak with secondary dominant intensity appears. Con- 
cerning this peak for h/t — 2.0, the values of energy w 
and residue z are {ui/t,z) = (3.95, 0.268), (3.96, 0.273), 
and (3.95, 0.265) for N = 16, 18, and 20, respectively. 
The weak size dependence of both values of to and z in- 
dicates that the secondary peak may be a delta-function 
contribution rather than a part of continuum spectrum 
in the bulk limit. 
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FIG. 4: Left: Pair spectral function with d a .2_ H 2-wave sym- 
metry in the 2D t-J-h model with 18 sites, J/t — 0.4, and 
various values of h/t. The delta functions (vertical bars) are 
broadened by a Lorentzian with a width of Q.lt (solid curves). 
Right: Spectral weight Z 2 h as a function of h/t in the 2D t-J-h 
model with J/t = 0.4. 



We estimate the spectral weight defined as 26 ' 27 
K* (7V h = 2)|A* ot |*o(^Vh-0))| 2 



%2h = 



(* (N h = 0)|(At ot )tAt ot |* (iV h = 0)) ' 



(4) 



which corresponds to the coherent peak of Pd(u>) at 
ui = 0. Note that Z 2 h is between and 1 because the 
denominator of Eq. (4) is equal to integration of Pd(^) 
over ui. The right panel of Fig. 4 shows i^h as a function 
of h/t. The weight is monotonically increasing function of 
h/t. Again, the size dependence of Z 2 h is weak for a large 
h/t. Therefore we expect that the coherent peak survives 
in the thermodynamic limit for a sufficiently large h/t. 

While our results suggest the enhancement of hole pair- 
ing by a staggered field, it is possible that the attraction 
between holes leads to phase separation. We calculate 
the clustering energy given by 20 

E c = E (N h = 4) + E (N h = 0) - 2E (N h - 2). (5) 

If this quantity is negative, the phase separation is ex- 
pected to occur. The results for N — 18 and 20 are shown 
in Fig. 5, which suggests that the region < h/t ^ 2 is 
not interrupted by the phase separation. 28 



The staggered field should induce a finite magnetic 
moment on each site, and consequently the commensu- 
rate AF long-range order. In Fig. 6 we show the result 
on the staggered-spin correlations given by C sp i n (r) = 
(WE; (-lY* +r HS?S? +r ) with r = (r x ,r y ). Indeed, 
the correlations seem to remain finite in the long-distance 
limit for h/t > 0, as expected. An important point is 



that both the d 2 .2_j / 2-wave superconducting correlations 
and the staggered-spin ones are enhanced by a staggered 
field. Therefore we expect the simultaneous presence of 
the d a; 2_j / 2-wave superconducting order and the commen- 
surate AF order in a certain range of h/t. 
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FIG. 5: Clustering energy as a function of h/t in the 2D t-J-h 




(I) 



t|oo i|t ^ tlaXolt — tfToolt 



(Ila) 

\ 1 1 1 1 



1 1 1 1 1 , 1 1 1 1 1 



t[oojjt 
1 1 1 1 1 

(lib) 

t \ 1 1 
t 



t aX-o t 



1 1 1 1 1 



1 1 1 1 



1 1 1 



t |oo| l 

1 1 1 1 



t 
\ 

1 1 1 1 



t[loo]t 
1 1 1 1 1 



tit* 
t 



I to 



t [To| l 
1 1 1 1 



FIG. 7: Typical processes of the second-order perturbation in 
the t-Jz-h model for h/t ^> 1. (I) is the ID case, while (Ha) 
and (lib) are the 2D case. The spin surrounded by a broken 
circle is against the direction of the staggered field. The three 
sites surrounded by a box correspond to (i,j,£) in the second 
term of Eq. (6). 



FIG. 6: Equal-time staggered-spin correlations as a function 
of distance in the 2D t-J-h model with N = 18, J/t = 0.4, 
and h/t = 0.0, 0.4, 1.0, and 2.0. The symbols are the same as 
in the left panel of Fig. 1. 



IV. PERTURBATION FROM LARGE-h/i LIMIT 

Why does the presence of a staggered field help the 
pairing formation? This may be understood from the 
l&vge-h/t limit. In order to treat analytically, we consider 
the t-J z -h model where the isotropic Heisenberg term (J) 
in Eq. (1) is replaced by the Ising term (J z ), following 
Refs. 14 and 15. For h/t > 1 we can regard the single- 
particle hopping term (t) as a perturbative one. In the 
Hilbert space with all spins along the direction of the 
staggered field, the second-order perturbation leads to 
the effective Hamiltonian given by Ti. c s = P7i. c sP where 



(i,i> 



-t 



-(other terms). 



Here and are the nearest neighbors, and 

p = nr GA ( 1 - ^YljeBi 1 - n ft)- The second term 
in the right-hand side of Eq. (6) includes hopping of a 
hole which occurs only if there is another hole in the 
neighboring site. Namely, it gives hopping of a hole pair. 
This is generated in the second-order processes shown in 
Fig. 7. In the 2D case the pair-hopping integral t is given 
by t 2 /(h + 3J z /2) where the denominator indicates the 
energy difference between the initial state and the inter- 
mediate one [see Fig. 7(IIa) and 7(IIb)]. We note that 
t is replaced by t 2 /(h + J z /2) in the ID case [see Fig. 
7(I)]. 14 < 15 



(6) 



The fact that the effective Hamiltonian includes the 
hole-pair hopping implies that hole binding and super- 
conductivity may be naturally realized in the antifer- 
romagnetic background forced by the staggered field. 
The effective pair- hopping integral t decreases as h/t in- 
creases. This suggests that a too large staggered field 
makes the hole pairs less mobile, which may be related 
to the phase separation at large h/t discussed above. On 
the other hand, it does not mean that a smaller h/t is 
better for superconductivity, because the calculation is 
based on the perturbation theory in t/h. 
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V. SUMMARY AND DISCUSSION 

We have investigated the binding energy and various 
correlation functions for the 2D t-J model in a staggered 
field. For the low-hole-density region at J/t = 0.4, the 
presence of a staggered field strengthens the attraction 
between two holes and helps the d x 2_ y 2-~wave supercon- 
ductivity. This implies that the commensurate antiferro- 
magnetic order and d a .2_ y 2-wave superconductivity can 
coexist in a strongly correlated electron system in two 
dimensions. 

One may ask whether calculation of small clusters with 
N ~ 20 provides some conclusive statements in a model. 
In fact, for the 2D t-J model without a staggered field, 
binding effects for N <~ 20 can be different from those for 
much larger size. 22,29 ' 30 However, the presence of a stag- 
gered field makes the coherence length (i.e., the size of a 
Cooper pair) small, and therefore the data for N <~ 20 is 
likely to reach the bulk limit for a sufficiently large stag- 
gered field, as evidenced by the weak size dependence. 
Thus our conclusion regarding the coexistence of super- 
conductivity and antiferromagnetism should hold at least 



near the boundary of the phase separation (h/t <~ 2). 31 
This would be of a conceptual interest, although such a 
large staggered field seems unrealistic. An open ques- 
tion is whether the picture for such a large field connect 
continuously to that for smaller field. For a realistic ap- 
plication of the present model, perhaps we need to know 
the effect of a small staggered field, as the effective stag- 
gered field produced by the interlayer coupling would be 
tiny. Unfortunately, for a small staggered field, the size 
dependence is still large and we cannot draw a definite 
conclusion from our present study based on small clus- 
ters. We hope future studies to clarify this question and 
its relation to the experiments. 
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